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ABSTRACT 

We estimate the metallicity in the broad emission-line region of the redshift z = 4.16 quasar, BR2248- 
1242, by comparing line ratios involving nitrogen to theoretical predictions. BR2248-1242 has unusually 
narrow emission lines with large equivalent widths, thus providing a rare opportunity to measure several 
line-ratio abundance diagnostics. The combined diagnostics indicate a metallicity of ~2 times solar. 
This result suggests that an episode of vigorous star formation occurred near BR2248-1242 prior to the 
observed z = 4.16 epoch. The time available for this enrichment episode is only ~1.5 Gyr at z — 4.16 
(for Ho = 65 km s _1 Mpc -1 , fi m = 0.3 and Q\ < 1). This evidence for high mctallicities and rapid star 
formation is consistent with the expected early-epoch evolution of dense galactic nuclei. 

Subject headings: galaxies: active — quasars: emission lines — quasars: individual 
(BR 2248-1242)— galaxies: formation 



1. INTRODUCTION 

The prominent emission line spectra of quasars provide 
valuable information on the physical state and chemical 
composition of the gas close to the quasars. If this gas was 
processed by stars in the surrounding host galaxy then the 
abundances we measure from the emission lines are a diag- 
nostic of the star formation history in these environments 
(see Hamann & Ferland 1999 for a general review). At 
redshifts z > 4, when the universe was only ~1 to 2 Gyr 
old (depending on cosmological models, see Figure 1 in 
Hamann & Ferland 1999), the results are especially inter- 
esting because they provide information about very early 
star formation, galaxy evolution, and chemical enrichment 
(Ostriker & Gnedin 1997; Friaca & Terlevich 1998). 

Ratios of emission lines involving nitrogen, N, are valu- 
able in determining the metallicity, Z, because of the 
expected "secondary" N production via the CNO cycle 
(Shields 1976). The nitrogen abundance therefore scales 
as Z 2 = (O/H) 2 and N/O scales as Z = O/H, providing 
a sensitive metallicity diagnostic even when direct mea- 
sures of Z = O/H are not available (Shields 1976; Ferland 
et al. 1996; Hamann & Ferland 1992, 1993, 1999). Ob- 
servations of H II regions indicate that secondary nitrogen 
production dominates for Z ;> 0.2 Z Q (Van Zee et al. 1998; 
Vila-Costas & Edmunds 1993). When secondary nitrogen 
production dominates, the N abundance is given by: 



[N/O] « [O/H] - q a log(Z/Z ) - q 



(1) 



where the square brackets indicate logarithmic ratios rel- 
ative to solar (Hamann et al. 2001). In environments 
with high star formation rates and short enrichment times, 
chemical evolution models predict that nitrogen produc- 
tion should be delayed relative to oxygen and carbon. In 



Equation 1 , q is the logarithmic offset caused by this delay. 
It is believed to range from ~0. 0-0.1 for "slow" chemical 
evolution, such as in H II regions, to ~0.2-0.5 for rapid 
evolution, such as may occur in massive galactic nuclei 
(Hamann et al. 2001). 

Studies of both emission and absorption lines in quasars 
suggest that their metallicities are near or above the solar 
value (Dietrich et al. 1999, 2000; Osmer et al. 1994; Con- 
stantin et al. 2001; Petitjean et al. 1994; Hamann 1997; 
Hamann & Ferland 1999; Hamann et al. 2001). These 
metallicities are consistent with both observational stud- 
ies (Pettini 2001) and theoretical simulations (Cen & Os- 
triker 1999) showing that metal abundances at any time in 
the universe are dependent on the local density. At every 
epoch, higher density regions, such as the central regions 
of galaxies, where quasars reside, should have much higher 
metallicities than lower density regions such as galactic 
halos or the intergalactic medium. This dependence of 
metallicity on density seems to be stronger than any rela- 
tionship to age (i.e. old does not necessarily mean metal 
poor). Evidence that density affects metallicity much more 
than age is found even in our own galaxy, where the central 
bulge contains many old but metal rich stars (Mc William 
& Rich 1994; Idiart et al. 1996). 

Much of the previous emission line analysis relied on 
the line ratios of N V A1240/He II A1640 and N V/C IV 
A1549. Ratios of various weak intercombination lines can 
be used to test that analysis. In this paper, we measure the 
emission lines in the rest frame UV spectrum of the high 
redshift (z = 4.16) quasar BR2248-1242 (Storrie-Lombardi 
et al. 1996), and we examine line ratios involving nitro- 
gen to estimate the metallicity. We selected BR2248-1242 
from a sample of 44 z J> 4 quasars in Constantin et al. 
(2001) because it has unusually narrow emission lines with 
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large equivalent widths, making the weak intercombina- ogy (http:/ /aeldata. phy.nist.gov/PhysR.efData/contents 



tion l jnes, such as N IV] A1486, O III] A1664, and N III] 
A1750, easier to detect. This is the most comprehensive 
emission-line abundance analysis so far for a z > 4 quasar. 

2. DATA & ANALYSIS 

The spectrum of BR2248-1242 was obtained at the 
MMT Observatory with the Red Spectrograph on Septem- 
ber 11, 1994. The full width at half maximum (FWHM) 
spectral resolution is —10 A which yields a difference of 
only —10% between the observed and actual FWHM of 
even the most narrow component. The wavelength range 
was chosen to span the redshifted Lya A1216 to N III] 
A1750 interval. This corresponds to an observed wave- 
length range of -5500 A to nearly 10,000 A. The observa- 
tion and data reduction details are discussed in Constantin 
et al. (2001). Figure 1 shows the observed spectrum. 

We measured the emission lines using tasks in the IRAF 5 
software package. We obtained the Galactic H I column 
density (Dickey & Lockman 1990) and calculated a redden- 
ing of Eb-v = 0.0693. We then performed a reddening 
correction and used the task NFIT1D to fit the continuum 
with a powerlaw of the form F v ex v a . The fit has a spec- 
tral index, a = —1.09 (see Figure 1) and is constrained 
by the flux in wavelength intervals between the emission 
lines, namely 6913-6991 A, 8733-8898 A, and 9278-9777 
A. We next used the task SPECFIT (Kriss 1994), which 
employs a \ 2 minimization routine, to fit each line with 
one or more Gaussian profiles. Figure 2 shows the fits. Ta- 
ble 1 lists the line fluxes and rest-frame equivalent widths 
(REWs) as measured above the fitted continuum. The 
fluxes, REWs, and FWHMs given in Table 1 are from the 
total fitted profiles, which can include several multiplet 
components and up to 3 Gaussian profiles per component 
(see below). Table 2 lists the rest wavelength, flux, REW, 
and FWHM of each individual Gaussian profile. We cal- 
culate the emission line redshift to be z = 4.156 based 
on the narrow Gaussian profiles fit to Lya and C IV. We 
use this redshift to calculate the rest wavelength for each 
component. 

We fit C IV A1549 with three Gaussian profiles, includ- 
ing a very broad profile (FWHM « 22500 kms" 1 ). This 
broad pedestal includes the unidentified —1600 A emission 
feature that has been noted in earlier studies (Laor et al. 
1994; Boyle 1990; Wilkes 1984). We ignore the doublet 
splitting in C IV because it is small compared to the ob- 
served line widths. N IV] A1486, He II A1640, and O III] 
A1664 are all fit with one Gaussian profile. Each com- 
ponent of the Si II AA1527, 1533 and Si IV AA1393, 1403 
doublets are fit with one Gaussian profile. The fluxes of 
these components are tied to a 1:1 ratio because this yields 
a better fit than a 2:1 ratio. In every multiplet fit, we fix 
the ratio of the central wavelengths to the known ratio 
of the rest wavelengths and tie together the FWHMs of 
each component. The N III] A1750 and O IV] A1403 mul- 
tiplets are fit with one Gaussian profile for each of their 
5 components. The flux component ratios in these cases 
are tied to the ratio of the statistical weight g times the 
A value for each transition, using atomic data obtained 
from the National Institute of Standards and Technol- 

5 The Image Reduction and Analysis Facility (IRAF) is distributed by 
the Association of Universities for Research in Astronomy, Inc. (USRA), 



atomic.html ) and Nussbaumer & Storey (1982). Lya 



A1216 is fit with three Gaussian components, with only 
the red side used for \ 2 minimization because of contami- 
nation by the Lya forest on the blue side. Each component 
of the N V A 1240 doublet is fit with two Gaussian profiles 
with FWHMs fixed to the values obtained from the nar- 
row and intermediate width profiles of our fit to C IV. 
This connection between C IV and N V can be justified 
because they are both high ionization lines with similar 
excitation and emission properties (Dietrich & Wilhelm- 
Erkens 2000). The relative fluxes of the two components 
of the N V doublet are set to 2:1 for each component, again 
because this yields a better fit than a 1:1 ratio. We fit the 
Si II AA1260, 1264, 1265 triplet with one Gaussian profile 
for each of the three components. Because the flux of Si II 
is so much smaller than Lya and N V, Si II is not well 
defined by the x 2 minimization. 

The primary uncertainty in our flux measurements is the 
continuum location. Many of the lines are blended, which 
introduces more uncertainty, particularly for N V in the 
wing of Lya. The broad component in the C IV fit intro- 
duces another uncertainty because it includes the — 1600 
A feature. We do not include this broad component in our 
calculated of flux ratios (§3 below) because it is believed 
to be an unrelated emission feature (Laor et al. 1994). We 
estimate the one a standard deviation of our measurement 
of the flux of Lya to be —10% based on repeated measure- 
ments with the continuum drawn at different levels. By 
the same method, we estimate the uncertainty in C IV, 
N V, He II, and O III] to be -15-20% and the uncertainty 
in the remaining lines with REW < 10 A to be —25%. 

3. RESULTS 

Table 3 lists various measured flux ratios and the metal- 
licities inferred from comparisons to the theoretical results 
from Hamann et al. (2001). Our preferred estimates of the 
metallicity are obtained from the model in Hamann et al. 
that uses a segmented powerlaw for the photoionizing con- 
tinuum shape (see Figure 3). This continuum shape is a 
good approximation to the average observed continuum 
in quasars (Zheng et al. 1997; Laor et al. 1997). The 
metallicity ranges in Table 3 represent the range of results 
obtained by comparisons to all three different continuum 
shapes calculated by Hamann et al. and shown in Figure 
3. These theoretical uncertainties are in addition to any 
uncertainty in the measured line strengths (§2). 

All of the line ratios in Table 3 yield a metallicity of 
Z k 1-3 Z Q . From a theoretical viewpoint, N III]/0 III] 
is the most reliable of the intercombination line ratios 
(Hamann et al. 2001). N V/He II is also a useful ra- 
tio, whereas N IV] /O III] and N IV] /C IV can be unre- 
liable because they are more sensitive to non-abundance 
effects and generally yield a wider range of results than 
N III]/0 III]. Based on which ratios are most accurately 
measured, and which are most reliable from a theoretical 
viewpoint (Hamann et al. 2001), we estimate the overall 
metallicity of BR2248-1242 to be roughly Z w 2 Z . 

It is important to note, however, that the actual metal- 
licity of BR2248-1242 may be 2-3 times higher than Z « 

the National Optical Astronomy Observatories, which is operated by 
under cooperative agreement with the National Science Foundation. 
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2 Zq because the theoretical models we use assume q « 
in Equation 1 . For environments with rapid chemical evo- 
lution, such as massive galactic nuclei, q « 0.2-0.5 may 
be more appropriate (Hamann et al. 2001). Therefore, 
the metallicity of BR2248-1242 is broadly consistent with 
previous emission line results of Z w 1-9 Zq derived for 
other quasars (Hamann ct al. 2001; Dietrich & Wilhelm- 
Erkens 2000; Hamann & Ferland 1999; Korista et al. 1998; 
Ferland et al. 1996). 

For any appropriate q value, BR2248-1242 may still 
yield a lower bound on the metallicity of luminous quasars 
because it appears to have an unusually low black hole 
mass (and so may reside in a lower mass, lower metallicity 
galaxy). We estimate the black hole mass in BR2248-1242 
to be M B h * 4.77 x 10 8 M Q , based on the measured 
FWHM of C IV and the formulae given in Kaspi et al. 
(2000) 



Rblr = (32.9±?;g) 
M = 1.464 x 10 5 



AL A (5100A) 
10 44 ergs s- 1 

^BLR 



0.700±0.033 



/ _Rblr_\ / 
\\t- days ) \10 3 km s 



"FWHM 



It — days 
(2) 

T ) 2 M © (3) 



where -Rblr is the radius of the broad line region for Hf3 1 
and -Rblr(C IV) « 0.5 Rblr{H(3) is the radius of the 
C IV region (Peterson 2001). The FWHM of C IV (1620 
km s _1 ) is represented by »fwhm and M is the mass of 
the black hole. We derive AL^(1450A) from the observed 
flux (Dietrich et al. 2002 in preparation) and use a pow- 
erlaw of the form F v oc v a with a = —0.4 to approximate 
the continuum shape and estimate ALa(5100A). We select 
a = —0.4 based on average quasar spectra from Vandcn 
Berk ct al. (2001) and Dietrich et al. (2002 in prepara- 
tion). The derived black hole mass in BR2248-1242 is un- 
usually low because of its unusually narrow emission lines. 
Previous studies suggest a correlation between black hole 
mass and the overall bulge/spheroidal component mass of 
the surrounding galaxy (Gebhardt et al. 2001; Ferrarese & 
Meritt 2000; Laor 2001; Wandel 1999). Together with the 
well known relationship between the mass and metallicity 
of galaxies (Faber 1973; Zaritsky et al. 1994; Jablonka et 
al. 1996), these results predict a relationship between the 
mass of the black hole and the metallicity surrounding the 
quasar (Warner et al. 2002 in preparation). BR2248-1242 
appears to have a low black hole mass, and so may be 
expected to have relatively low metal abundances. 

4. DISCUSSION 

The result for Z £ Z in BR2248-1242 suggests that an 
episode of rapid and extensive star formation occurred be- 
fore the observed z = 4.16 epoch. The spectral similarity 
of BR2248-1242 to other z > 4 quasars (apart from the 
narrow line widths, see Constantin ct al. 2001) suggests 
further that the resulti for Z ;> Z Q should apply gener- 
ally. The time available for the chemical enrichment at 
z = 4.16 is only ~1.5 Gyr (for H = 65 km s _1 Mpc -1 , 
TO = 0.3 and f^A <, 1). This star formation episode early 
in the history of the universe would be consistent with 
observations showing that the central regions of massive 
galaxies today contain old, metal rich stars (McWilliam 
& Rich 1994; Bruzual et al. 1997; Worthey et al. 1992; 
Idiart et al. 1996). This high rcdshift evolution would 



also be well within the parameters derived in some recent 
simulations, which show that the densest proto-galactic 
condensations can form stars and reach solar or higher 
metallicities at z ;> 6 (Ostriker & Gnedin 1997; Haiman & 
Loeb, 2001). Recent observations of the highest redshift 
quasars suggest that reionization occurred around z ~ 6 
(Becker et al. 2001; Djorgovski et al. 2001). If the reion- 
ization is due to stars, this would provide further evidence 
for substantial star formation beginning at z ;> 6. Quasar 
abundance studies can provide observational tests of these 
evolution scenarios. The next step is clearly to extend the 
abundance analysis to the highest possible redshifts, and 
to compare the results across wide ranges in redshift, host 
galaxy types, and central black hole mass. 
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Fig. 1 — 
extinction. 



Spectrum of BR2248-1242 with our fit to the continuum overplotted (dashed line). The spectrum has been corrected for 




Fig. 2. — Gaussian fits of emission lines of BR2248-1242. The continuum is normalized by the power law fit shown in Figure f . The 
continuum (at unity) and composite line fits are plotted ag solid curves T>m individual f"i"r"'"'"' c of the fits are shown as dotted and dashed 
lines. Larger versions of these plots are available at ittp://www. astro. ufl.edu/~warner/BR2248/ 
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Fig. 3. — Comparison of flux ratios measured in BR2248-1242 to theoretical results from Figure 5 of Hamann ct al. (2001). Comparisons 
of our measured line ratios to the theoretical model represented by the solid line yield our preferred estimates for metallicity. The solid 
line is obtained from the model in Hamann et al. that uses a segmented powerlaw for the photoionizing continuum shape, which is a good 
approximation to the average observed continuum in quasars. The dotted line corresponds to the model that uses an incident spectrum 
defined by Matthews & Ferland (1987), while the dashed line corresponds to a powerlaw with index a = —1.0 across the infrared through 
X-rays . 



Table 1 
Emission Line Data 



Line 


A obs 


Flux 6 


Flux/Lya 


REW a 


FWHM C 


Lya A1216 


6270.3 


6040 


1.000 


246 


1450 


N V A1240 


6396.6 


420 


0.070 


18 


2380 


Si II A1263 


6509.9 


90 


0.014 


4 


3660 


I A1303 


6728.7 


250 


0.041 


10 


3020 


C II A1335 


6874.2 


30 


0.005 


1 


1270 


Si IV A1397 


7204.5 


180 


0.029 


8 


3490 


IV] A1403 


7254.8 


140 


0.023 


7 


3240 


N IV] A1486 


7663.9 


130 


0.021 


6 


1910 


Si II A1531 


7892.1 


120 


0.020 


6 


2450 


C IV A1549 


7984.8 


2400 


0.398 


122 


1620 


He II A1640 


8455.3 


210 


0.035 


11 


2200 


III] A1665 


8580.6 


350 


0.059 


19 


4200 


N III] A1751 


9029.2 


160 


0.026 


9 


3530 



In units of A 

In units of 10 -17 ergs cm~ 2 s 
In units of km s _1 



Table 2 
Emission Line Components 



Component 


^■rest 


Flux fc 


REW a 


FWHM C 


Lya 1 


1216.1 


1580 


64 


970 


Lya 2 


1218.1 


2780 


112 


2810 


Lya 3 


1225.1 


1690 


69 


10020 


N V 1 


1240.2 


140 


6 


1280 


N V 2 


1240.6 


280 


12 


3710 


Si II 


1262.8 


90 


4 


3390 


OI 


1305.0 


250 


10 


3020 


CII 


1333.2 


30 


1 


1230 


Si IV 


1397.3 


180 


8 


1570 


OIV] 


1407.1 


140 


7 


3000 


NIV] 


1486.3 


130 


6 


1910 


Si II 


1530.7 


120 


6 


1140 


CIV 1 


1548.6 


1070 


55 


1280 


C IV 2 


1552.2 


1330 


68 


3710 


C IV 3 


1551.0 


1820 


93 


22570 


He II 


1639.9 


210 


11 


2200 


OIII] 


1664.2 


350 


19 


4200 


N III] 


1751.2 


160 


9 


3390 



In units of A 

In units of 10 -17 ergs cm~ 2 s 
In units of km s _1 



Table 3 
Emission Line Ratios 



Lines 


Ratio 


z/z 


N III]/0 III] 


0.45 


i 74+0-16 


N IV] /O III] 


0.36 


9 ,,0+0.00 


N IV] /C IV 


0.05 


1 cc+0.16 
1 - oo -0.44 


NV/He II 


2.02 


o n7 +3.94 
•J-"' -1.48 


NV/C IV 


0.18 


1 40+0 09 



